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H , energies between 500 and 1500 GeV. We consider rather light Higgs scalars, with 

mass below the W^W^ threshold, produced in the double Higgs-strahlung channel 
e~^e~ — > HHZ and decaying via the dominant mode H — > hb. Assuming ideal vertex 
tagging of the heavy-quark jets and mass reconstruction of the Z boson produced in 
association with the Higgs pair, we compare the yield of the above process to those 
of all irreducible electroweak and strong backgrounds proceeding through a bbbbZ 
intermediate state. Total cross sections and differential spectra of phenomenological 
relevance to the selection and analysis of the signal are given and discussed. 



Abstract 

We assess the possibility of measuring the trilinear self-coupling of the Standard 
Model Higgs boson at future electron-positron linear colliders with centre-of-mass 
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1 Introduction 

In all probability a Higgs boson will be discovered at either LEP II, Run II of the Teva- 
tron, or the Large Hadron Collider (LHC). It is then inevitable that the emphasis of Higgs 
physics will be turned away from discovery and instead will focus on the investigation of 
Higgs boson properties, such as its mass, width and branching ratios. Although much 
interesting Higgs phenomenology can be done at the LHC, many analyses are made infea- 
sible by the rather messy nature of hadron colliders. Instead one must resort to the much 
cleaner environment of e^e~ annihilations, for example, at the Linear Collider (LC) [jl|, 
where precision measurements at the TeV scale can be made. 

One particularly interesting task to be carried out at future colliders is the recon- 
struction of the Higgs potential itself, possibly confirming, or denying, the mechanism of 
spontaneous electroweak symmetry breaking. This can be achieved by measuring the tri- 
linear Xhhh and quadrilinear Xhhhh Higgs self-couplings, which can then be compared 
with the predictions of the Standard Model (SM), or indeed the Minimal Supersymmetric 
Standard Model (MSSM)[[ 

A measurement of the trilinear term, Xhhh, is the first step in reconstructing the Higgs 
potential. At a future e^e^ collider, the Xhhh coupling of the SM is accessible through 
double Higgs-strahlung off Z bosons, in the process e"*'e~ —^ HHZ. This is the mechanism 
with which we will be concerned in this paper (for the MSSM see Ref. |§|, 0). The SM 
signal of the trilinear Higgs self-coupling has been thoroughly investigated in Ref. |^ (with 
its MSSM counterparts), and was found to be small but measurable for an intermediate 
mass Higgs boson, given a high integrated luminosity. In contrast, the quadrilinear vertex, 
Xhhhh, is unmeasurable at the energy scale of the proposed LCs due to suppression by 
an additional power of the electromagnetic coupling constant. 

However, in measuring Xhhh, one must be sure that the already small signal can 
be distinguished from its backgrounds without being appreciably reduced. Here we will 
examine the H —>■ bb decay channel over the Higgs mass range Mh ^ 140 GeV and present 
kinematic cuts to aid its selection. The case of off- and on-shell H — > H^ =*=(*) VT^ decays 
for Mh ^140 GeV is under examination elsewhere . 



■^In principle, the former coupling is amenable to investigation also at hadron and high energy photon 
colliders too, via double Higgs-strahlung off W^ or Z bosons H, D, W~^W~ or ZZ fusion [||, @, ||, 
gluon-gluon fusion Q or 77 fusion |^, ^, |^. 



If one assumes very efficient tagging and high-purity samphng of b quarks, the back- 
grounds to a Xhhh measurement from double Higgs events in the 46 decay channel are 
primarily the 'irreducible' ones via bbbbZ intermediate states |jlll, which can be separated 



into EW and QCD backgrounds. Furthermore, the double Higgs-strahlung process (see 
Fig. D: 

e+e" -^ HHZ -^ bbbbZ (1) 

contains diagrams proceeding via an HHZ intermediate state but not dependent on Xhhh 
(graphs 1-3 in Fig. 0), as well as that sensitive to the trilinear Higgs self-coupling (graph 
4 in Fig. |I|). In addition, we also include four extra diagrams, which differ only in the 
exchange of the four-momenta and helicities of two identical b quarks (or, equivalently, 
antiquarks) and a minus sign (due to Fermi-Dirac statistics pertinent to identical fermionic 
particles). However, the narrow width of the Higgs resonance ensures that the interference 
will be negligible and these extra diagrams could be included by symmetry. 

The other two backgrounds proceed via purely EW interactions (see Figs. §-§), 

e+e" -^ EW graphs -^ bbbbZ, (2) 

and via QCD couplings as well (see Fig. ^), 

e+e" -^ QCD graphs -^ bbbbZ, (3) 

and both contain no more than one intermediate Higgs boson. The EW background, 
process (0), is of 0{al^) away from resonances, but can, in principle, be problematic due 
to the presence of both Z vectors and Higgs scalars yielding bb pairs. Finally, the QCD 
background, process (^, is of 0{al^al) away from resonances. Here, although there are no 
heavy objects decaying to bb pairs, the production rate itself could give difficulties due to 
the presence of the strong coupling. As with process (|l|), one must include diagrams with 
the interchange of the two identical b (anti)quarks also in the EW and QCD background 
processes. In contrast to the signal, here interference effects are sizable. 

The plan of the paper is as follows. The next Section details the procedure adopted in 
computing the relevant scattering amplitudes. Sect. |^ displays our numerical results and 
contains our discussion. Finally, in the last Section, we summarize and conclude. 



2 The matrix elements (MEs) 

The double Higgs-strahlung process ([l|) proceeds at lowest-order through the diagrams of 
Fig. |l|, as explained in the Introduction. They are rather straightforward to calculate in 
the case of on-shell HHZ production (see, e.g., Ref. for an analytic expression of the 
ME). 

The EW background (Q) derives from many graphs: 550 in total (again, considering 
the b (anti)quark statistics). However, they can conveniently be grouped into different 
'topologies': that is, collections of diagrams with identical (non-)resonant structure. We 
have isolated 23 of these, and displayed them in Figs. |^ and 0, depending on whether one 
or zero Higgs intermediate states are involved, respectively. There are 214 graphs of the 
first kind and 336 of the second. This approach, of splitting the ME in (no n-) resonant 
subprocesses, facilitates the integration over the phase space and further provides an insight 
into the fundamental dynamics. On the one hand, one can compute each of the topologies 
separately, with the appropriate mapping of variables, thus optimizing the accuracy of the 
numerical integration. On the other hand, one is able to assess the relative weight of the 
various subprocesses into the full scattering amplitude, by comparing the various integrals 
with each other. However, one should recall that the amplitudes squared associated to 
each of these topologies are in general non-gauge invariant. In fact, the latter is recovered 
only when the various (non)-resonant terms are summed up. For reasons of space, we will 
not dwell in technicalities any further here, as a good guide to this technique can be found 



in Ref. |T^ . (The resonant structure of the various subchannels ought to be self-evident 
in Figs. ^.) 

The QCD diagrams associated to process (^ can be found in Fig. ^ In total, one has 
120 of these, with only five different (non-) resonant topologies. The integration in this 
case is much simpler than in the EW case and can in fact be carried out with percent 
accuracy directly over the full ME using standard multichannel Monte Carlo methods. 

Non-zero interference effects exist between processes (0), (0) and (^). However, given 
the very narrow width of the Higgs boson (always below 20 MeV over the mass range 
considered here), any interference with the signal can be safely neglected. Furthermore, 
we will see that the dominant subprocesses of the two backgrounds have very different 
topologies, so one also expects their interference to be negligible. Therefore, given that 
their calculation would be rather cumbersome, we do not consider them in our analysis. 



The large number of amplitudes can easily and efficiently be dealt with in the numerical 
evaluation if one resorts to helicity amplitudes. In doing so, we have adopted the HELAS 



subroutines [1^. The algorithm used to perform the multi-dimensional integrations was 



VEGAS 14 



Numerical inputs were as follows. The strong coupling constant Ug entering the QCD 
process (0) has been evaluated at two loops, with Nf = 5 and Ajj^ =160 MeV, at a scale 
equal to the collider CM energy, i/i = Ecm- The EM coupling constant was Ogm = 1/128. 
The sine squared of the Weinberg angle was sin^ 9w = 0.232. The fermionic (pole) masses 
were rrie = and m^ = 4.25 GeV. As for the gauge boson masses (and their widths), 
we have used: Mz = 91.19 GeV, Tz = 2.50 GeV, Mw = MzcosOw ~ 80 GeV and 
Tw = 2.08 GeV. 

Concerning the Higgs boson, we have spanned its mass Mh over the range 100 to 
150 GeV and we have computed its width, Th, by means of the program described in 



Ref. [0, which uses a running b mass in evaluating the H ^ hh decay fraction. Thus, 
for consistency, we have evolved the value of mi, entering the Hhh Yukawa coupling of the 
H —>■ bb decay currents in the same way. 

We have adopted as CM energies typical for the LC the values Ecm = 500, 1000 and 
1500 GeV. 

Notice that, in the remainder of this paper, total and differential rates are those at the 
partonic level, as we identify jets with the partons from which they originate. In order 
to resolve the latter as four separate systems, we impose the following acceptance cuts: 
E{b) > 10 GeV on the energy of each b (anti)quark and cos(6, b) < 0.95 on the relative 
separation of all possible 26 combinations. We further assume that b jets are distinguishable 
from light-quark and gluon jets (e.g., by using yU- vertex tagging techniques). However, no 
efficiency to tag the four b quarks is included in our results. Also, the Z boson is treated as 
on-shell and no branching ratio (BR) is applied to quantify its possible decays. In practise, 
in order to simplify the treatment of the final state, one may assume that the Z boson 
decays leptonically (i.e., Z — > i^i~, with i = e,fj,,T) or hadronically into light quark jets 
(i.e., Z -^ qq, with q j^ b). 

Finally, we have not included Initial State Radiation (ISR) [l^ in our calculations. In 
fact, we would expect it to affect rather similarly the various processes (|lD-(0). As we are 
basically interested in relative rates among the latter, we are confident that the salient 
features of our results are indifferent to the presence or not of photons radiated by the 



incoming electron-positron beani^. 



3 Results 



The total cross sections for process (|T]), at the three CM energies considered here, can be 
found in the top-left frame of Fig. |^, as a function of Mh- The decrease of the total rates 
with increasing Higgs mass is mainly the effect of the BR of the decay channel H — * 66, see. 



e.g.. Fig. 1 of Ref. ||T5[. This mode is dominant and very close to 1 up to the opening of the 
off-shell H — i> 14^='=*14^'t: decay channel, which occurs at Mh ~ 140 GeV. In contrast, the 
production cross section for e^e~ -^ HHZ is much less sensitive to Mh 0]. In addition, 
because reaction (|1|) is an annihilation process proportional to 1/s, a larger CM energy 
tends to deplete the production rates, as long as Ecm ^ 2Mh + Mz- When this is no 
longer true, e.g., at 500 GeV and Mh ^ 140 GeV, phase space suppression can overturn 
the 1/s propagator effects. This is evidenced by the crossing of the curves for 500 and 
1000 GeV. 

In practice, the maximum cross section for double Higgs-strahlung (|^) is reached at 
energies E^m ~ 2Mh + Mz + 200 GeV 0. For Higgs masses in the lower part of the 
Mh range considered here, e.g., Mh = 110 GeV (where the bottom-antibottom channel 
is unrivaled by any other decay mode), this corresponds to Ecm = 500 GeV. Furthermore, 
the sensitivity of the production rates of reaction (|T]) on Xhhh is higher at lower collider 
energies Q. Thus, in order to illustrate the interplay between reactions (|l])-(§), we will 
in the following focus on the case of a CM energy of 500 GeV, top-right corner of Fig. ^, 
as an illustrative example. In fact, the discussion for the other two choices, E^m = 1000 
and 1500 GeV (the two bottom plots of Fig. ^), would be rather similar, so we refrain 
from repeating it. (Also note that the signal-to-background {S/B) ratio improves with 
increasing energy.) 

The rise at 500 GeV of the purely EW background (0) with the Higgs mass can be 
understood in the following terms. The dominant components of the EW process are those 
given by: 

1. e'^e^ -^ ZZZ -^ bbbbZ, first from the left in the second row of topologies in Fig. ^. 
That is, triple Z production with no Higgs boson involved. 



^We also neglect beamsstrahlung and Linac energy spread, by assuming a narrow beam design M&. 



2. e"'"e^ -^ HZZ -^ bbbbZ, first(first) from the left(right) in the fifth(fourth) row of 
topologies in Fig. ^ (also including the diagrams in which the on-shell Z is connected 
to the electron-positron line). That is, single Higgs-strahlung production in associa- 
tion with an additional Z, with the Higgs decaying to bb. The cross sections of these 
two channels are obviously identical. 

3. e~^e~ -^ HZ -^ Z*Z*Z -^ bbbbZ, first from the right in the third row of topologies 
in Fig. ^. That is, single Higgs-strahlung production with the Higgs decaying to bbbb 
via two off-shell Z* bosons. 

4. e~^e^ -^ ZH — » bbZ*Z -^ bbbbZ, first(first) from the right(left) in the first(second) 
row of topologies in Fig. |^. That is, two single Higgs-strahlung production channels 
with the Higgs decaying to bbZ via one off-shell Z* boson. Also the cross sections of 
these two channels are identical to each other, as in 2. 

The production rates of 1.-4. as separate subprocesses can be found in the upper portion 
of Fig. ^. All other EW subprocesses are much smaller and rarely exceed 10~^ femtobarns, 
so we do not plot them here. 

The QCD process (^ is dominated by e^e^ -^ ZZ production with one of the two Z 
bosons decaying hadronically into four b jets. This subprocess corresponds to the topology 
in the middle of the first row of diagrams in Fig. ^. Notice that Higgs diagrams are involved 
in this process as well (bottom- right topology in the above figure). These correspond to 
single Higgs-strahlung production with the Higgs scalar subsequently decaying into bbbb via 
an off-shell gluon. Their contribution is not negligible, owning to the large ZH production 
rates, as can be seen in the lower portion of Fig. ^. The somewhat unexpected dependence 
of the latter upon Mh (with a maximum at 130 GeV) is the result of the interplay between 
our acceptance cuts and phase space effects. The contribution of the other diagrams, which 
do not resonate, is typically one order of magnitude smaller than the ZZ and ZH mediated 
graphs, with the interferences even smaller (and generally negative). 

One should note from Fig. ^ that the overall rates of the signal are quite small (also 
recall that we neglect tagging efficiency as well as the Z decay rates), even at low Higgs 
masses where both the production and decay rates are largest. In fact, they are always 
below 0.2 femtobarns for all energies from 500 to 1500 GeV, although this can be dou- 
bled simply by polarizing the incoming electron and positron beams. Thus, as already 
recognised in Ref. 0, where on-shell production studies of process (p were performed. 



luminosities of the order of one inverse attobarn need to be collected before statistically 
significant measurements of Xhhh can be performed. This emphasizes the need of high 
luminosity at any future LC 

We now proceed by looking at several differential spectra of reactions (|l])-@, in order 
to find suitable kinematic cuts which will enhance the S/B ratio. The distributions in 
E{b) and cos (6, b) leave little to exploit in separating signal from background after the 
acceptance cuts are made, especially with respect to the EW background. We turn then 
to other spectra, for example, invariant masses of b (anti)quark systems. In this respect, 
we have plotted those of the following combinations: 

(a,b) of 26 systems, for the case in which the b jets come from the same production vertex 
('right' pairing) and the opposite case as well ('wrong' pairing); 

(c) of 36 systems, in which only two 6 (anti) quarks have the same EM charge; 

(d) of the 46 system. 

We denote the mass spectra of the systems (a)-(d) as Mr (66) and M^ibb) (where R{W) 
signifies the right (wrong) combination), M{bbb) and M{bbbb), respectively. In the first 
three cases, there exists more than one combination of 6 quarks. In such instances, we bin 
them all in the same distribution each with identical event weight. Further notice that the 
26 invariant masses that can be reconstructed experimentally are actually appropriately 
weighted superpositions of M^(66) and Mw{bb). In particular, if the 6 charge tag is 
available, then it is roughly the sum of the two. If not, the latter is about twice as large 
as the former. 

The invariant mass spectra can be found in Fig. ^, for the combination Ecm = 500 
GeV and M^ = 110 GeV. Here, one can appreciate the narrow Higgs peak^ in the Mji{bb) 
distribution, that can certainly be exploited in the signal selection, especially against the 
QCD background, which is rather fiat in the vicinity of Mh- Apparently, this is no longer 
true for the EW process, as it also displays a resonance at Mh (induced by the diagrams 
in Fig. 1^ which carry an external on-shell current H ^ bb). However, events of the type 
(ffl) contain two 26 invariant masses naturally peaking at Mh, whereas only one would 



^Recall that for AIh = 110 GeV one has Th ~ 3 MeV. The Higgs resonances in the top-left frame of 
Fig. have been smeared out by incorporating a 5 GeV bin width, emulating the finite efficiency of the 
detectors in determining energies and angles. 



appear in samples produced by process (0) (apart from accidental resonant mispairings) . 
Thus, even in the case of the EW background one can achieve a significant noise reduction. 
Finally, requiring that none of the 26 invariant masses reproduce a Z boson will also be 
helpful in this respect, as evident from the Mji{bb) spectrum of process (0). However, in 
this case, the invariant mass resolution of di-jet systems must be at least as good as the 
difference [Mfj — Mz)/2, in order to resolve the Z and H peaks. Other mass distributions 
can be useful too in reducing the noise while keeping a substantial portion of the signal. 
Of some help are the M{bbb) and M{bbbb) spectra. In particular, notice that the minimum 
value of the latter is about 2Mh for process (|l|), whereas for reaction (0) it is lower, 
typically around 2Mz or Mh + Mz, as driven by the two dominant components of the 
EW background at low Higgs mass (i.e., subprocesses 1. and 2., respectively, see top of 
Fig. ^. The QCD background can stretch to M{bbbb) values even further below the 2Mh 
end point of the signal (the more the larger Mh), indeed showing a peak both at Mh and 
Mz, corresponding to the H ^ 4b and Z —^ Ab decays induced by the second and last 
topologies in Fig. ^. As for the M{bbb) spectrum, its shape is strongly related to that of 
the 46 mass. In a sense, by excluding one of the four b quarks from the mass reconstruction 
corresponds to smearing the M{bbbb) distributions, so that the broad prominent peak at 
M{bbb) ^ 90 — 100 GeV in the case of the QCD process can be viewed as the superposition 
of the remains of the two narrow ones seen in M{bbbb). For this very reason then, once a 
selection cut is imposed on one of the two masses, this is very likely to affect the other in 
a similar manner. 

The spectacular differences seen in Fig. |^ between, on the one hand, process (j^), and, 
on the other hand, reactions (0) and (^ (more in the former than in the latter), are a 
direct reflection of the rather different resonant structure of the various channels, that is, 
the form of the time-like propagators (i.e., s-channels) in the corresponding MEs. However, 
one should expect further kinematic differences, driven by the presence in the backgrounds 
of space-like propagators (i.e., t,-u-channels), which are instead absent in the signal (see 
Figs. |1]-§|). This is most evident in the properties of the four-quark hadronic system (d) 
recoiling against the Z boson. As the internal dynamics of the four b quarks is very 
different in each process (|l|)-(|l), we also study the cases (a)-(c) separately. 

One can appreciate the propagator effects by plotting, for example, the cosine of the 
polar angle (i.e., with the beam axes) of the four b quark system (or, indeed, the real 
Z). See Fig. |, where, again, E^m = 500 GeV and Mh = 110 GeV. Notice that the 
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backgrounds are much more forward peaked than the signal. This can be understood by 
recalhng that the QCD events are mainly due to e~^e~ -^ ZZ production followed by 
the decay of one of the gauge bosons into four b quarks. The ZZ pair is produced via 
t, M-channel graphs, so the gauge bosons are preferentially directed forward and backward 
into the detector. In contrast, the signal (|l|) is always induced by s-channel graphs. The 
EW background (Q) has a more complicated structure but is still sizably dominated by 
forward production. The behaviour of cos{bbb) and cos{bb) is very similar to that of the 
four-quark system. In practise, it is the strong boost of the Z bosons produced forwards 
and backwards in processes (^-(0), combined with the small value of rrif, (compared to 
the typical process scale, i.e., -Ecm), that produces a similar angular pattern for all multi-6 
systems of the background, regardless of their actual number. This is true for reaction (^ 
also. Therefore, all angular distributions displayed can boast strong (though correlated) 
discriminatory powers, allowing one to separate signal and backgrounds events efficiently. 

An alternative possible means of disentangling the effects of the propagators is to resort 
to the differential distributions of the above systems (a)-(d) in transverse momentum, pT. 
These are plotted in Fig. P, for the same E^m and M^ as the previous two figures. However, 
this kinematic variable proves not to be useful. In fact, the only discriminating distribution 
is the pT for all four b quarks (or equivalently the final state Z boson), and this only singles 
out the QCD background, a large fraction of which populates the range beyond 180-200 
GeV. Neither the signal nor the EW background do so and always look rather similar in 
their shape (even in the spectra involving 26 quarks only, once these are appropriately 
combined together). 

Notice that two of the b quarks in the QCD background originate from gluon splitting 
and are therefore rather soft. Consequently, one expects the four b quarks to be more 
planar, in the 46 rest frame, for the QCD background than for the signal, where they are 
all the decay products of heavy bosons. We study this by plotting the thrust T fl^ and 



sphericity S [|1^ distributions in Fig. 0. Indeed, such quantities could prove useful in 



reducing the QCD background but are harmless against the EW background. 

Before proceeding to apply dedicated selection cuts, we remark that kinematic features 
similar to those displayed in Figs. |7HlO| can also be seen at the other two values of CM 



energies considered here and for other Higgs masses. In fact, increasing the value of 
Mff in those distributions mainly translates into a 'movable' resonant peak in Mji{bb) as 
well as lower-end point in M{bbbb) and into somewhat softer (harder) spectra in invariant 



mass (transverse momentum) than at the smaller Mh value considered so far. Moreover, 
the t, M-channel dependence of the backgrounds, as opposed to the s-one of the signal, is 
more marked at higher Ecm values. Finally, for angular distributions, a larger Higgs mass 
does not remove the big differences seen between the three channels ([l|)-(0). 

Therefore, in all generality, following our discussions of Figs. f^-ITOl, and recalling the 
need to economize on the loss of signal because of its rather small production and decay 
rates, one can optimise the S/ B ratio by imposing the cuts: 

\M{bh) — Mh\ < 5 GeV (on exactly two combinations of 26 systems), 

\M{bb) — Mz\ > 5 GeV (for all combinations of 26 systems) 

M{bbbb) > 2Mh, \ cos(26, 3b, 46) | < 0.75. (4) 

In enforcing these constraints, we assume no b jet charge determination. Moreover, the 
reader should recall that the spectra of the four hadronic systems (a)-(d) are all corre- 
lated, in each of the quantities studied above, and so are the invariant masses, transverse 
momenta and polar angles among themselves. 

The counterpart of Fig. ^ after the implementation of the above cuts is Fig. |rT[ The 
effect of the latter is a drastic reduction of both background rates (@)-(^, while maintain- 
ing a large portion of the original signal (p. Further notice how the imposition of the cuts 
(^) modifies the hierarchy of cross sections for process (|I|) with the CM energy, as now the 
largest rates occur at Ecm = 1000 GeV and the smallest at Ecm = 500 GeV (compare to 
Fig. D. 

The S/B ratios turn out to be enormously large for not too heavy Higgs masses. 
For example, at E^^ = 500(1000) [1500] GeV and for Mh = 110 GeV, one gets S/B = 
25(60) [104], where S corresponds to the rates of reaction (|^) and B refers to the sum of 
the cross sections for processes (|2|)-(^. The reduction of both backgrounds amounts to 
several order of magnitudes, particularly for the case of the strong process, whereas the 
loss of signal is much more contained. The acceptance of the latter is better at higher 
collider energies and lower Higgs masses. In fact, the poorest rate occurs for Ecm = 500 
GeV at the upper end of the Mh range, where more than 90% of the signal is sacrificed. 
We should however remark that the suppression of the backgrounds comes largely from the 
invariant mass cuts on M^f, advocated in (^. (In fact, they are crucial not only in selecting 
the Mh resonance of the signal, but also in minimizing the rejection of the latter around 
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Mz when inispairings occur: notice the shoulder at 90 GeV of the Mw{hh) spectrum of 
reaction (|T])). The value we have adopted for the resolution is rather high, considering the 
large uncertainties normally associated with the experimental determination of jet angles 
and energies, though not unrealistic in view of the most recent studies [1^ . The ability of 



the actual detectors in guaranteeing the performances foreseen at present is thus crucial 
for the feasibility of dedicate studies of double Higgs-strahlung events at the LC. 

A related aspect is the efficiency of tagging the h quarks necessarily present in the final 
state of reaction (||), particularly in the case in which the Z boson decays hadronically. 
On the one hand, given the high production rate of six jet events from QCD pO[ and 



multiple gauge boson resonances |21| in light quark and gluon jets, it is desirable to resort 
to heavy fiavour identification in hadronic samples. On the other hand, the poor statistics 
of the HHZ signal requires a judicious approach in order not to deplete the latter below 
detection level. According to recent studies [^, the two instances can be combined 
successfully, as efficiencies for tagging hb pairs produced in Higgs decays were computed 
to be as large as e^^i ~ 90%, with mis-identification probabilities of light (charmed) quarks 
as low as eqq(^cc) ~ 0.3(4)% (and negligible for gluons). If such a projection for the LC 
detectors proves to be true, then even the requirement of tagging exactly four h quarks 
in double-Higgs events of the type (||) might be statistically feasible, thus suppressing the 
reducible backgrounds to really marginal levels [^ . 

One must also bear in mind that experimental considerations, such as the performances 
of detectors, the fragmentation/hadronization dynamics and a realistic treatment of the 
Z boson decays, are also important when determining what cuts should be made. Such 



considerations are beyond the scope of this paper, and are under study elsewhere ||Tl |. 

Finally, the number of signal and backgrounds events seen per inverse attobarn of 
luminosity at -Ecm = 500, 1000, and 1500 GeV, with Mh = 110 GeV, can be seen in 
Tab. m. Of course, one could relax one or more of the constraints we have adopted to try 
to improve the signal rates without letting the backgrounds become unmanageably large. 
For example, by removing the cuts on cos{bb) and cos{bbb) one can enhance the signal 
rates by about a factor of two. However, the EW background would also increase by a 
comparable amount and the QCD rate would do by a somewhat larger factor, of at least 



three/four. Kinematic fits can also help in improving the S/B ratios [11] 
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4 Summary 

In conclusion, the overwhelming irreducible background from EW and QCD processes of 
the type e~^e~ -^ bbbbZ to double Higgs production in association with Z bosons and decay 
in the channel H -^ bb, i.e., e'^e" -^ HHZ -^ bbbbZ, should easily be suppressed down to 
manageable levels by simple kinematics cuts: e.g., in invariant masses and polar angles. 

The number of signal events is generally rather low, but will be observable at the LC 
given the following 'mandatory conditions' (some of which have already been outlined in 



Ref. ITTI): 

• very high luminosity; 

• excellent b tagging performances; 

• high di-jet resolution. 



The requirement advocated in Ref. |TT| of a good forward acceptance for jets may also be 



added to the above list, as we have explicitly verified (though not shown) that single jet 
directions can stretch in the double Higgs-strahlung process up to about 20 degrees in polar 
angle. Finally, beam polarization can also be invoked to increase the signal-to-background 
rates 0. 
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Number of Events per ab ^ after selection cuts 


Ecm = 500 GeV 


Ecm = 1000 GeV 


Ecm = 1500 GeV 


signal 


26 


40 


34 


Electroweak 


1.0 


0.6 


0.3 


QCD 


0.032 


0.026 


0.016 



Table 1: The number of events for signal and backgrounds per inverse attobarn of lumi- 
nosity after selection cuts for centre-of-mass energies of 500, 1000 and 1500 GeV, a Higgs 
mass of 110 GeV, and with polarized electron and positron beams. 
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Figure 1: Diagrams contributing at lowest-order to e^'e^ -^ h:ibJ)^hQZ^ via purely EW 
interactions containing two Higgs bosons in intermediate states. An internal wavy line 
corresponds to a Z boson (labeled as Z). The total number of actual diagrams is 4. Finally, 
diagrams which differ from those above only in the exchange 3 ^-* 5 (or, equivalently, 4 ^^ 6) 
must also be considered, preceded by a minus sign. 
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Figure 2: Topologies contributing at lowest-order to e^e^ ^ h'^^^^Z-j via purely EW 
interactions containing only one Higgs boson in intermediate states. An internal wavy line 
corresponds to a 7 or a Z (labeled as A and Z, respectively), as appropriate. The total 
number of actual diagrams is 107. Finally, diagrams which differ from those above only 
in the exchange 3 ^^ 5 (or, equivalently, 4 ^^ 6) must also be considered, preceded by a 
minus sign. 
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Figure 3: Topologies contributing at lowest-order to e^'e^ -^ h^4b^QZ-j via purely EW 
interactions containing no Higgs bosons in intermediate states. An internal wavy line 
corresponds to a 7 or a Z (labeled as A and Z, respectively), as appropriate. The total 
number of actual diagrams is 168. Finally, diagrams which differ from those above only 
in the exchange 3 ^^ 5 (or, equivalently, 4 ^^ 6) must also be considered, preceded by a 
minus sign. 
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Figure 4: Topologies contributing at lowest-order to e^'e^ -^ h^hih^Jj^Z'j via QCD interac- 
tions containing one gluon in intermediate states. An internal wavy line corresponds to 
a 7 or a Z (labeled as A and Z, respectively), as appropriate, whereas a helical one refers 
to a g. The total number of actual diagrams is 60. Finally, diagrams which differ from 
those above only in the exchange 3 ^^ 5 (or, equivalently, 4 ^^ 6) must also be considered, 
preceded by a minus sign. 
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Figure 5: Top-left: cross sections in femtobarns for the signal at three different collider 
energies: 500, 1000 and 1500 GeV. Top-right (Bottom-left) [Bottom-right]: cross sections in 
femtobarns for the signal versus the EW and QCD backgrounds at 500(1000) [1500] GeV. 
Our acceptance cuts in energy and separation of the four b quarks have been implemented. 
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Figure 6: Top: cross sections in femtobarns for the four dominant components (see the 
text) of the purely EW background. Bottom: cross sections in femtobarns for the total 
and Higgs components (see the text) of the QCD background. The CM energy is 500 GeV. 
Our acceptance cuts in energy and separation of the four h quarks have been implemented. 
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Figure 7: Differential distributions in invariant mass of multi-jet systems containing one, 
two, three and four h (anti) quarks. The CM energy is 500 GeV and the Higgs mass 
110 GeV. Our acceptance cuts in energy and separation of the four h quarks have been 
implemented. 
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Figure 8: Differential distributions in (cosine of the) polar angle of multi-jet systems 
containing one, two, three and four h (anti)quarks. (The spectrum of the 46 system is 
identical to that of the Z boson.) The CM energy is 500 GeV and the Higgs mass 110 GeV. 
Our acceptance cuts in energy and separation of the four h quarks have been implemented. 
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Figure 9: Differential distributions in transverse momentum of multi-jet systems containing 
one, two, three and four h (anti)quarks. (The spectrum of the 46 system is identical to 
that of the Z boson.) The CM energy is 500 GeV and the Higgs mass 110 GeV. Our 
acceptance cuts in energy and separation of the four h quarks have been implemented. 



^ 
-^ 



b 



e e ^ X ^ 4b Z 
E_ = 500 GeV 



solid ^ HHZ, Mh = HO GeV 
dashed ^ EW 



dotted ^ QCD 
(no interferences) 



T 1 1 1 1 1 1 1 r 



T 1 1 r 



J I I ^ 




0.5 



0.6 



0.7 



0.8 



0.9 



M 



m 
b 



3.0 

2.5 

2.0 

1.5 

1.0 K 

0.5 



T 1 1 1 1 1 1 1 r 



T 1 1 r 



, I A / I 




M 

fD 



^ 







0.2 



0.4 



0.6 



0.8 



Figure 10: Differential distributions in thrust and sphericity. The CM energy is 500 GeV 
and the Higgs mass 110 GeV. Our acceptance cuts in energy and separation of the four b 
quarks have been implemented. 
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Figure 11: Top-left: cross sections in femtobarns for the signal at three different collider 
energies: 500, 1000 and 1500 GeV. Top-right (Bottom-left) [Bottom-right]: cross sections in 
femtobarns for the signal versus the EW and QCD backgrounds at 500(1000) [1500] GeV. 
Our acceptance cuts in energy and separation of the four b quarks have been implemented 
along with the selection cuts of eq. (H). 



